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Abstract 

Recently a new picture has been developed for examining Wilson lines, and the corresponding 
anomalous dimensions which govern their renormalization properties. By making a particular 
coordinate transform, the calculation of the cusp anomalous dimension in QED or QCD can be 
related to the energy of a pair of static charges in Euclidean Anti-de-Sitter (AdS) space. This 
paper shows how the same picture can be used to describe Wilson lines in quantum gravity. 
We show how the relevant cusp anomalous dimension (which has recently been shown to be 
one loop exact) can be obtained using the Newtonian limit of General Relativity. We also show 
how both the QED and gravity cases emerge as special cases of a general formulation, and that 
a continuous parameter exists which interpolates between them. The results may be useful in 
examining the relations between gauge and gravity theories. 



1 Introduction 

Wilson lines have been studied for many years in a variety of contexts [THlO|. in both Abelian 
and non-Abelian gauge theories. In particular, they govern the structure of large logarithms in 
perturbation theory due to soft gluon emission |llU12j as is also seen in other approaches such as 
soft collinear effective theory (SCET) [13-20J, factorisation theorems |21H23| . or the path integral 
technique of [24J. Essentially, hard partons emitting soft gluons cannot recoil and thus can only 
change by a phase. For this phase to have the right gauge transformation properties to slot into an 
amplitude, it must be a Wilson line evaluated along a contour given by the hard momentum of the 
outgoing particle. The ultraviolet renormalisation properties of Wilson lines govern the infrared 
singularities of scattering amplitudes [5T-f7t l25H28] , a fact which has recently been used to derive an 
all-order ansatz for the infrared singularities of QCD scattering amplitudes [29-33J, whose conse- 
quences have been discussed further in |34H37j . More exotic applications of Wilson lines come from 
the interface between gauge and string theory, where it has recently been conjectured that certain 
Wilson loop configurations are dual to scattering amplitudes in M = 4 Super- Yang Mills theory at 
strong (and weak) coupling [3H|l39]. 
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Correlators of Wilson line operators have ultraviolet (UV) divergences. For smooth contours, these 
are completely removed by renormalisation of the gauge theory coupling. For contours which con- 
tain cusps or intersections, however, additional UV singularities are present which, after renormal- 
isation of the coupling, factor out into an overall multiplicative factor [3H2] . Wilson line correlators 
are then governed by renormalisation group equations, involving an anomalous dimension. For the 
case of a single cusp, the latter is known as the cusp anomalous dimension, which has become a 
quantity of central importance in the study of infrared singularities and resummation. 

In Abelian theories, the all-order structure of the cusp anomalous dimension is dictated by a simple 
subset of possible Feynman diagrams, namely those in which the Wilson line contours meeting at 
the cusp are joined by a connected subgraph. The only possibilities are single photon emissions, 
or photon emissions which are joined by fermion loops. Traditionally this result was derived using 
the eikonal identity for multiple soft photon emissions [21J. Recently, an alternative proof has 
been given, using path integral techniques that relate the exponentiation of connected subgraphs 
in Wilson loops with a cusp, to the textbook exponentiation of connected diagrams in a general 
quantum field theory [23]. This simple structure is referred to as Abelian exponentiation in the 
literature, and in particular implies that if there are no propagating fermions, the cusp anomalous 
dimension is one-loop exact. Consequently, the infrared singularities of scattering amplitudes are 
purely governed by the exponentiation of the one- loop result. A similar result holds for correlators 
involving more than two Wilson lines meeting at a point. Then, the relevant anomalous dimension 
is completely determined to all orders by connected subdiagrams spanning the outgoing particles. 
Again, one-loop exactness follows if fermion loops are absent. 

In non-Abelian theories, the cusp anomalous dimension is, unsurprisingly, more complicated. Nev- 
ertheless, its all order structure is still dictated by a subset of Feynman diagrams, which permit a 
simple topological classification: as well as connected subdiagrams spanning the outgoing Wilson 
line contours, one may also have diagrams which are two-eikonal line irreducibl^. Such diagrams 
have modified colour factors which are maximally non-Abelian, and are referred to as webs in the 
literature [3QH320. Another way to think about webs is that they are, by definition, those di- 
agrams which enter the exponent of the Wilson line correlator. Their structure is thus usually 
referred to as non-Abelian exponentiation, by analogy with the Abelian case. These results were 
recently reconsidered in [24], where an alternative derivation of webs was given using statistical 
physics methods. Furthermore, the notion of webs has recently been extended to correlators in- 
volving many Wilson lines meeting at a point, by two groups of authors [44,45]. By analogy with 
the two line case (alternatively, a single contour with a cusp), webs are those diagrams which enter 
the exponent of the Wilson line correlator. However, their structure is markedly different to the 
two line case. The most significant difference is that reducible diagrams contribute to the exponent 
in general. Furthermore, webs become closed sets of diagrams whose members are related by gluon 
permutations. Their colour and kinematic information is entangled to all orders in perturbation 
theory by so-called web mixing matrices [45], whose structure is in principle governed purely from 
combinatorics [46]. This structure has yet to be fully explored, but is certainly crucial to under- 
standing the all-order structure of infrared singularities in multileg scattering amplitudes [47] . 



3 That is, subdiagrams which cannot be disconnected by a single cut through both eikonal lines. 

4 A novel geometric interpretation of webs in two parton scattering has very recently been presented in [43] . 
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Recently, a new way to think about Wilson lines was presented in |48| . The authors start off by 
considering a collection of Wilson lines meeting at a point in Minkowski space. They then transform 
to a set of radial coordinates, which map Minkowski space to M x AdS, where AdS denotes Eu- 
clidean anti-de-Sitter space. Parametrisation invariance of the Wilson line contours in Minkowski 
space becomes time translation invariance in the transformed space, such that each Wilson line can 
be considered as a static charge. In addition, the dilatation operator in Minkowski space (whose 
eigenvalues are dimensions) maps to the operator d T in radial coordinates, where r is the time-like 
coordinate normal to the (space- like) slices of AdS. Eigenvalues of the latter operator are energies 
(up to a factor of i), so that one may obtain the anomalous dimension associated with a given 
Wilson line correlator by evaluating the energy of a collection of static charges in Euclidean AdS 
space. This is worked out in detail in [38] for QED and QCD at one loop, and the conceptual 
picture thus obtained is used to motivate a family of conformal gauges for use in field theory Wil- 
son line calculations. At two loops, as shown explicitly in [48] . such gauges eliminate non-Abelian 
graphs containing a three-gluon vertex. Although such graphs are not completely absent at higher 
orders, it seems likely that such gauges may provide a natural framework for probing the conjec- 
tured all-order structure of IR singularities in QCD fixed- angle scattering amplitudes [30, 321133]. 
which involves correlations between only pairs of partons in the case that they are all massless. 

The purpose of this paper is to extend the radial picture for Wilson lines to perturbative quan- 
tum gravity, in the form of General Relativity (GR) minimally coupled to (scalar) matter. The 
consideration of infrared singularities in gravity dates back to [49] . and there has recently been 
renewed interest [501152] . based on a number of motivations. Firstly, although GR contains non- 
renormalisable ultraviolet singularities, alternative field theories of gravity may have the same long 
distance behaviour [53] . Secondly, there is a growing body of work involving intriguing connections 
between scattering amplitudes in gauge and gravity theories, at arbitrary loop level (see e.g. [54H56] . 
and [57U58] for recent applications), which may have a string theoretic origin. There then exists 
the real possibility that knowledge about quarks and gluons can tell us about gravity, or vice versa. 
To this end, it is important to develop common conceptual pictures, that allow us interpret similar 
physics (such as infrared singularity structures) in a natural way, an observation which motivates 
the present study. 

The structure of the paper is as follows. In section [2] we review the radial picture for Wilson 
lines developed in [48] . as well as recent results from perturbative quantum gravity, which will be 
necessary for the rest of the paper. In section [3j we apply the radial picture to Wilson line operators 
in quantum gravity, showing explicitly that the form of the one-loop cusp anomalous dimension in 
gravity can be reproduced from the Newtonian limit of General Relativity in Euclidean AdS space. 
In section [4j we comment on the lightlike limit of gravitational Wilson line operators, in particular 
observing the cancellation of collinear singularities in the radial picture. In section [5] we introduce 
a general formulation of the cusp anomalous dimension calculation, and show that a continuous 
parameter exists that interpolates smoothly between the QED and gravity cases. In section [6l we 
discuss our results before concluding. Some technical details are collected in appendices. 
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2 Review of necessary concepts 
2.1 The radial picture for Wilson lines 

In this section, we briefly review the approach of [H] for describing Wilson lines as static charges 
in AdS space. Given that we will be focussing explicitly on gravity in the rest of the paper, we here 
discuss only QED, thus ignoring additional complications due to non-trivial colour structure. 



Let us start with the definition of a Wilson line operator in an Abelian gauge theory: 



W(C) = exp 



ie I dx ll A^(x) 



(1) 



where A^{x) is the gauge field, e the coupling constant, and the line integral is taken over the 
contour C in Minkowski space. In applications of Wilson lines to scattering amplitudes, one is 
interested in Wilson line contours which are fully determined by the momenta of the outgoing 
particles. Each Wilson line then represents an outgoing particle dressed by an infinite number 
of soft photon emissions, and the soft part of the scattering amplitude is given by the vacuum 
expectation value of the product 



r poo 

W(ni, . . . , n L ) = TT exp ie \ dstpi ■ A(spi) 

i=i L Jo 



(2) 



where L is the number of hard outgoing particles. Here we have parameterised each (straight) 
Wilson line contour by Xi = SiPi, where pi is the 4-momentum of outgoing particle i. Note that 
then has dimensions of (length) 2 . It is more common in the literature to instead write 
where n« is the 4-velocity. The parameter tj then has dimensions of length, but one is free to make 
the above choice, given that the Wilson line operator is invariant under rescalings 



As,; 



(3) 



The fact that we choose to use the momentum rather than the velocity in eq. (pQ) is to make ex- 
plicit the analogy between gauge theory and gravity, where in the latter theory momenta are more 
important given that they play the role of charges. 



Consider now Minkowski space with time coordinate t, and spherical polar coordinates r, 9 and 
4>, representing the radial, polar and azimuthal coordinates respectively. Then a given Wilson line 
direction can be parametrised by [28] 

= e T (cosh p, sinh (3n) , (4) 

where n is a unit 3- vector. In these so-called radial coordinates, the Minkowski space metric 
becomes^] 

ds 2 = e 2r [dr 2 - (dji 2 + sinh 2 /3 d£l\)\ , (5) 

where dd 2 . is the squared line element on a 2-sphere. Next, one interprets — oo < r < oo as a time 
coordinate, and uses the fact that Abelian gauge theory (in the absence of propagating fermions) is 

5 Note that in this subsection only we use the metric (+,-,-,-), for ease of comparison with [48]. Throughout the 
rest of the paper, we will use the alternative choice (-,+,+,+), which is more common in the gravity literature. 
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classically conformally invariant in four dimensions. That is, one is free to rescale the line element 
of eq. © by an overall factor, and thus to replace eq. © by 



ds 2 = dr 2 - (dp 2 + sinh 2 pdSli). (6) 

The time coordinate is now explicitly decoupled from the spatial coordinates, and spatial slices 
(i.e. at fixed r) constitute Euclidean AdS space in three dimensions. Two important identifications 
between the original Minkowski space and the transformed space are as follows: 

• Reparametrisations of the Wilson line contour in Minkowski space, e.g. the rescalings of 
eq. ([3]), map to time translations in the transformed space. Thus, reparametrisation invariance 
of Wilson lines shows up, in the latter space, as invariance under shifts in r. This allows one 
to consider the Wilson lines in the transformed space as static charges. 

• The dilatation operator x^d^ (whose eigenvalues are dimensions) maps to the operator d r in 
the transformed space. This is related to the Hamiltonian operator in the transformed space 
via 

d T = iU RxAdS , (7) 

so that the anomalous dimension of a Wilson line correlator maps to the energy of the cor- 
responding collection of static charges in the Euclidean AdS space, with an additional factor 
of i. 

The electrostatic potential </> due to the two charges is given by the solution of Laplace's equation 
in Euclidean AdS space: 

V ^ = -Arz a f> ( sinh2 = °> ( 8 ) 



sinh" p 



whose solution is 



cf>{P) = Ci + C 2 coth/3, (9) 

where the {Ci} are constants of integration. As explained in detail in [48J, this does not have 
the correct behaviour as ft — > oo, where the cusp anomalous dimension should diverge linearly in 
P, due to the appearance of collinear singularities in Minkowski space. This is due to the fact 
that the solution implicitly includes the effect of a spurious charge corresponding to a phantom 
initial state particle (i.e. eq. ([9]) has a pole at P = and P = in). One may correct for this 
by first analytically continuing the spacelike part of the radial coordinate space to a Euclidean 
3-sphere, and subsequently adding a constant charge density to the source term of the equation for 
the electrostatic potential. These constant densities cancel out upon constructing any collection of 
Wilson line operators, which must satisfy charge conservation. The corrected solution is 

<t>= 4^2 [(7r + i/3)coth/3 + C], (10) 

where the overall normalisation is fixed by the amount of constant charge necessary to cancel the 
phantom charges. 

Having constructed the potential, the total energy of a pair of charges (which in [48] is obtained 
by first constructing the electric field) is 

E{p 12 ) = f^r [(*■ + ifa) coth/3i 2 + C] , (11) 
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where P12 is the geodesic distance between the charges, which is identified with the cusp angle in 
Minkowski space. The constant C can be fixed by analytic continuation to the situation in which 
one of the Wilson lines is incoming. The vanishing of the cusp anomalous dimension when the 
two Wilson lines become parallel then implies that C = —i [48]. Finally, one finds that the cusp 
anomalous dimension should be given by 

r = iEifa) = [(*7T - £12) coth p 12 + 1] , (12) 

in exact agreement with an explicit Minkowski space field theory calculation. The above argument 
is for the case of QED. The generalisation to a non-Abelian context is straightforward, and is also 
presented in [48J. 



2.2 Infrared singularities in Quantum GR 

In the previous section, we have reviewed the arguments of [48J which relate the cusp anomalous 
dimension in a gauge theory to the energy of a static charge configuration in Euclidean AdS space. 
The aim of the present paper is to outline how the same arguments can be applied in quantum 
gravity. To this end, we must first recap salient features regarding the structure of infrared singu- 
larities in GR. For ease of comparison with previous literature on quantum gravity, we will from 
now on adopt the (-,+,+,+) metric, in contrast with the previous subsection. 

Here we will consider general relativity minimally coupled to a scalar field <3>, a theory whose total 
action is given by 

S = S E . H .[g^} + S m ,t[^^,9 ll % (13) 
where Se.h. ^ s ^ ne Einstein-Hilbert action, and 

S mat [<*>*, <n = J d d x^~g [-grdp&h* - m 2 <T<I>] (14) 

in d dimensions, where g is the determinant of the metric tensor g^ u . Perturbation theory can be 
defined by expanding the metric tensor according to 

9nu = V(j,v + Kh ^u, (15) 

where hnv is then the graviton field, and k = y/WirG^, with Gm Newton's constant. Note that 
a different choice is often made in the literature (see e.g. [59]), in which one instead defines the 
graviton via the expansion of the quantity 

sT = yf=9<T- (16) 

This simplifies the Feynman rules for emission of gravitons from scalar particles, and was the choice 
adopted by [51] which derived Wilson line operators for soft graviton emission. Here, however, we 
will stick with the choice of eq. (|15p . for reasons that will become clear. Substituting eq. (|15p into 
eq. p4p and expanding up to first order in k (including also the factor y/—g), the matter action 
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becomes^ 



S mat [$*,<S>,hn = j d d x 



-(■;"<i> ^,<i> - ™ 2 $*$ + k ( -^V"*** - -7/"^$*^$ 



+ 0(k 2 ). (17) 



The momentum-space Feynman rule for single graviton emission from a scalar is then 

j [(-m 2 +Pl ■p2W+P f iP2+PiPZ] , (18) 

where p\ and P2 are both outgoing. We will also need the graviton propagator, which in d dimensions 
is given (in the de Donder gauge) by [60] 



Given the theory defined by eq. (|13p . one may consider scattering amplitudes involving a num- 
ber L of external scalars, which interact by the exchange of gravitons. When emitted gravitons 
become soft, infrared divergences occur. The study of infrared singularities in quantum general 
relativity dates back to the classic paper |49j, which established that Abelian exponentiation can 
be generalised to gravity from the QED case. Earlier this year, attention again turned to the IR 
sector of gravity, with a particular emphasis on describing the singularity structure using methods 
developed in the context of gauge theory [50H52] , leading both to an extension of our knowledge of 
IR effects in gravity, and also a more unified view of the shared physics underlying the IR sectors 
of gravity and gauge theories. It has been known for some years that amplitudes in gauge theory 
with L external lines factorise into the schematic form [6TT-I65I 



A L = s-n- JJJi. (20) 



i=i 



Here S is a soft function, collecting all infrared singularities due to soft gauge boson emission, Ti 
is a hard function which is finite in dimensional regularisation as d — > 4, and Jj is a jet function 
collecting hard collinear singularities associated with the i th external lin^Zl- The soft function is 
given by a vacuum expectation value of Wilson lines meeting at a point, where there is one Wilson 
line for each external particle. The authors of [50] suggested that a similar structure should hold 
for pure quantum gravity, namely that one has 

A L = S ■ H (21) 

in the gravitational context, where again S and H denote soft and hard functions, but in this 
case there are no jet functions, due to the fact that collinear singularities cancel after summing 
over diagrams. The latter fact was first established in [49], for soft-collinear emissions. As in the 



6 Note that after expanding in the weak field approximation, contractions involving upper and lower indices are 
interpreted as involving the Minkowski metric ?7 M „. 

7 Strictly speaking, there is a double counting in eq. (|20p. given that soft-collinear singularities are present in both 
the jet functions and the soft function. This can be easily rectified by dividing by eikonal jet functions J7i, which we 
may ignore in the present discussion. 
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Figure 1: The emission of a graviton from an external scalar line. 



gauge theory case, the soft function is given by a Wilson line correlator, where an appropriate 
form of the Wilson line operator (representing the emission of soft gravitons from an eikonal line) 
was suggested in [500. Furthermore, the anomalous dimension governing the renormalization of 
gravitational Wilson line correlators was argued to be one loop exact. Reference [51] investigated 
some of these issues further in the context of GR minimally coupled to matter (as described by 
eq. (I13p ). and derived explicit forms for the Wilson line operators of quantum gravity using path 
integral techniques borrowed from gauge theory [23]. This also generalised the discussion to include 
the effects of massive eikonal lines, and corrections beyond the eikonal approximation. Finally, the 
authors of [52] have undertaken a complete and rigorous analysis of the infrared singularity struc- 
ture of pure gravity amplitudes, using techniques earlier applied to QCD (see [66J for a pedagogical 
review). This confirmed the one- loop exactness of the anomalous dimension underlying the IR 
singularity structure of gravity amplitudes. In addition, the cancellation of collinear singularities 
was generalised from soft to hard collinear configurations. 

For this paper, we will need the form of the gravitational Wilson line operator for the weak field 
expansion of eq. (|15p . for a straight line contour along the direction of an outgoing particle. This 



for momentum p^. Note that this is not the same as the result presented in [51], owing to the 
fact that that paper used a weak field expansion based on eq. (|16p . However, we can justify this 
result as follows. First, note that the emission of a graviton from an external scalar line is given by 
figure [TJ This introduces a propagator for the intermediate scalar line, and an emission vertex for 
the graviton. The latter is given by the vertex of eq. (|18p with p\ = —(p — k) and P2 = P (recalling 
that the {p^} were defined to be outgoing above), which in the soft limit k — > becomes 



is given by [H 




(22) 



- [-(p 2 + m 2 )rf v - 2p^p v ] 



= —inp^p u . 



(23) 



Combining this with the propagator for the intermediate scalar line gives 




(24) 



8 Note that what we refer to as a Wilson line in this paper is not to be confused with the parallel transport operator 
of GR, defined in terms of the Christoffel symbol, which is also sometimes referred to in Wilson line terms. 
9 An alternative Wilson line operator was studied in [67]. See also [50] for a discussion of this point. 



The right-hand side constitutes an effective Feynman rule for the emission of a soft graviton from an 
eikonal line. In particular, it is insensitive to the spin of the emitting particle. This same Feynman 
rule is indeed generated by the Wilson line operator of eq. (|22p . which shows that the latter is the 
correct operator. 

An interesting consistency check of the above Wilson line operator is that one may derive Newton's 
law of gravity, by calculating the expectation value of a certain Wilson loop in Minkowski space. 
This is analogous to the derivation of the Coulomb potential in QED. Although the derivation of 
Newton's law from perturbative quantum gravity is by no means a new result [60J, this calculation 
is usually performed by taking the non-relativistic limit of the full one-graviton exchange graph. 
Given that the alternative derivation gives extra insight into the gravitational Wilson line opera- 
tors, we present this here in appendix [AfPl . 

Armed with the above results, we are now ready to examine Wilson line operators in quantum 
gravity, using the Euclidean AdS picture of [48]. This is the subject of the following section. 

3 Gravitational Wilson lines in the radial picture 

In the previous sections, we have reviewed both the radial coordinate picture for calculating cusp 
anomalous dimensions of Wilson line operators, and also the known structure of infrared divergences 
in quantum GR. In this section, we combine these concepts, in order to elucidate the gravitational 
analogue of the analysis of [H] for abelian and non-abelian gauge theories. As stated in the previous 
section, we use the metric (-,+,+,+), so that the equivalent of the radial coordinate space of eq. ([5]) 
is 



Our starting point is to note that, unlike the case of conventional gauge theories, the Wilson line 
operator of eq. (|22p is not invariant under rescalings of the form of eq. ([3]). This does not prevent 
us from examining Wilson lines in radial coordinates. However, it is not then true that the system 
in the transformed space will be invariant under time translations (in r). In particular, one must 
include the prefactor of e 2r in the metric of eq. (|25|) . which changes the nature of computations 
involving energies in radial space. In calculating the latter, we will be interested in the component 
h TT of the graviton field. The lack of rescaling invariance in Minkowski space implies that this will 
be r-dependent. We can surmise this dependence by rewriting the Wilson line operator of eq. (I22p 
as 



where symmetry of the graviton field (hav = h Ufl ) means that we can decide to take either of the 
momentum factors out of the integral along the Wilson line parameter. The exponent now has the 
explicit form of a reparametrisation-invariant term 



ds 2 = e 2r [-dT 2 + {dp 2 + sinh 2 /3dO|)]. 



(25) 




(26) 




(27) 



This calculation is closely related to a similar analysis in [B5], as we explain in the appendix. 



9 



multiplied by a charge factor p^. That this charge is the 4-momentum of the emitting particle is 
as expected from general relativity. Rewriting the expression (|27p as 



oo 







drx v h„ v (x^, (28) 



where the dot represents differentiation with respect to r, it follows that 



= Yt Kt = °' (29) 



so that the invariance under r translations results from the invariance of the expression (|28p under 
reparametrisations. That is, the quantity h TV is independent of r in the radial coordinate space 0. 
As a consequence, one may write 

h TT = x^x v h liu = e T K(p), (30) 

where we have used the fact that x^ = x^ = e T n^, with independent of r. Given that is in 
the direction of the hard momentum p^ of the particle emitting soft gravitons, we may write 

n» = !-, x» = e T — , (31) 
m m 

where the normalisation ensures n 2 = —1. 

Analogously to the QED analysis of [48J, one may consider calculating the potential energy of a 
pair of charges in the radial coordinate space. In this case, these will be masses, which we label 
by mi and mi- The total potential energy can be constructed by first considering the potential 
energy generated by the mass located at the origin of the radial coordinate space, and then 
considering the mass mi as a test particle a geodesic distance /3 away. 

To this end, we must construct the equation satisfied by h TT in radial coordinate space, and the 
first step is to note that the Wilson line phase of eq. (|26|) can be written as 

in „ r°° ik ro ° 



-p v \ dx» h^(s^) = -p v drh TU . (32) 

1 JO 1 J~oo 

Inserting a delta-function in the radial coordinate space, this can be further rewritten as 



-P 



/OO P 
dr / d 3 x5 {3) (x)h TU . (33) 
-OO J 



2 

Equating this to the general form of a source term in curved space, 

d A x^f u K 



one finds that the current which sources the conformally invariant quantity h TV is given by 

^r = ^p v 5^{x). (34) 



11 Care is needed here in interpreting the notation. The index v in h TV is a Minkowski space index, which has yet 
to be transformed to the radial coordinate space. 
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One then finds 

Ur=9rrjJ l = -^p^\ X ), (35) 

where we have absorbed geometric factors into the delta function so that this is now normalised 
according to 

r <fxJ^8^{x) = 1, (36) 



where is the determinant of the spatial part of the metric g^ u , disregarding the overall factor 
of e 2r . The component h TT of the graviton field satisfies the wave equation (see e.g. [69J) 

\Jh TT =j TT , (37) 

where from eqs. (|3T1 1351) one finds 

3rr = ^j^ T = ^m 2 e- T 5^{x). (38) 
The left-hand side of eq. (I37p is the covariant D'Alambertian operator 

□ = (V=^U) . (39) 

We see from eq. (|37p that gravitational Wilson lines in the radial space correspond to time- 
dependent sources. This is a manifestation of the overall e 2r factor in the metric of eq. f)25[) : 
the radial size of the spacelike part of the space changes with r. If we consider a mass fixed at 
the origin, the flux of the gravitational field through an expanding surface of radius ~ e T (3 will be 
constant. If we instead consider a surface at fixed (5 around the origin, the flux through such a 
surface grows as a function of time, such that it looks as if a mass me~ T is enclosed* 2 !, which is 
precisely the content of eq. (|3j 



Using the metric of eq. (|25l) . eq. ([371) becomes 

e -2r ( V 2 _ 3) Kr = ^ m2e - 5 (3) {x) t ( 40) 

where V 2 is the three-dimensional LaplaciarJ^l. The homogeneous equation is not simply Laplace's 
equation as in the case of QED, but rather the Helmholtz equation. The additional term linear in 
h TT is a direct consequence of the fact that the flux through a surface of fixed /3 is changing with 
time due to the overall time dependence of the metric, as discussed above. The relation of the 
graviton field to the Newtonian potential is (again see e.g. [69] ) 

$ = \Kt> (41) 

so that eq. ()40j) may be rewritten as 

V 2 $ - 3$ = ^ m 2 e T 5 {3) (x). (42) 

12 One must also remember that there is a factor \f—g involved in calculating the flux, so that a mass which 
decreases with time is consistent with a flux which grows with time. 

13 Note that, for ease of comparison with the QED case, we denote by V 2 the Laplacian without the r dependence 
implied by the metric of eq. (|25ll . 
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This is the gravitational analogue of Laplace's equation for the QED case of [48J. It is essentially 
Newton's law of gravity, but with an extra term proportional to <!> resulting from the use of a 
time-dependent curvilinear coordinate system. Prom eqs. {|30|, |4*T|) we may write 

(43) 



$(P,t) = K([3)e\ K(J3) = -Ktf). 



Then using the Laplacian of eq. (JSj> , we must find the general solution of the homogeneous equation 



which is 



sinh /3 



K(P) = A, 



dp sinh^ PdpK(/3) - 3K((3) = 0, 



1 



+ 2 sinh /3 



+ A 2 cosh p. 



(44) 



(45) 



sinh (5 

Note that, as in the QED case, the general solution of the homogeneous equation is a superposition 
of a function which is even under the transformation of eq. (150p . and a function which is odd. These 
are the first and second terms in eq. (|45p respectively. As shown in appendix [Bj one can fix A\ by 
considering a 3-surface consisting of a cylinder of radius /3q and height (in time) To- The result is 

,^2 



K(f3) = m 2 



16vr 



1 



sinh j3 



+ 2 sinh f3 



+ A 2 cosh/3. 



(46) 



Another feature shared with the QED case is that this solution has the wrong behaviour as f3 — > oo. 
As shown in appendix £3 the actual result ~ fie" at large /3. This is a factor e@ up on the QED 
case (where the divergence is linear), due to the fact that there are extra momentum factors in 
gravity. The failure of the solution (|46p to reproduce this behaviour is, as in the QED case, due 
to a spurious charge. Noting that eq. (|46p has a pole at both f3 = and j3 = in, this solution 
represents the combined effect of the physical mass at (3 = 0, and a spurious mass at f} = in. 

In [48] , and as reviewed in section (2[ this problem is solved in the QED case by modifying the 
current j T by a constant charge density. A similar procedure can be used in the gravity case, as 
we now discuss. First, recall that the current is given by eq. (I35h . and that this sources the 
conformally invariant quantity involving one Minkowski-space and one radial-space index. If 
we were to solve for hn T , we would find a charge at the origin (3 = of the radial space, and 



ITT. 



a spurious charge — (corresponding to an incoming momentum in Minkowski space) at /3 
The analogue of adding a constant charge density in the present case is to modify the current of 
eq. (f35l) so as to give 



JfJ.7 



,-2t 



(47) 



where K is such as to remove the spurious incoming momentum. Any given collection of gravita- 
tional Wilson lines will obey momentum conservation, such that the constant charge densities thus 
added will cancel out. This is the analogue of the cancellation of the constant terms due to electric 
charge conservation in the QED case. 



Above, we solved directly for h TT , which is related to the Newtonian potential in the radial coordi- 
nate space. To this end, one must consider the modified current 



Jt 



K 

2 6 ~ 



5 {?) {x) + K 



(48) 



12 



Here = x^ = e T (cosh /3, sinh/3n) is a general point in the radial coordinate space. In the first 
term, this is constrained to be related to the momentum (corresponding the test charge at the 
origin) by eq. ([3T]) . due to the delta function. In the second term this is not the case, and using 



771-2(1,0) (in Minkowski coordinates) one finds 



K 

r 



m 2 



5 {3) (x) + K cosh f3 



(49) 



That is, the effect of a constant current density for the conformally invariant quantity hu T becomes 
a non-trivial charge density ~ cosh (3 in the current for h TT . That this makes physical sense can be 
seen as follows. We have already noted that the spurious pole of eq. (f46l) corresponds to a mass at 
(3 = itt, in addition to the mass at /3 = 0. Both of these masses are necessarily positive. However, 
the charge density cosh f3 has the property of being odd under the transformation 



iir 



(50) 



Thus, the added charge density reinforces the physical charge in the upper branch of the AdS space, 
but acts to cancel out the fake charge in the lower branch, as required. Absorbing various factors 
in the constant charge density K, eq. (|44p is modified to 



1 



sinh (3 



dp sinh 2 pdpK(P)) - 3K((3) = B cosh (3, 



(51) 



whose general solution is 



K(J3) = (A 1 + A 3 f3) 



sinh (3 



+ 2 sinh (3 



+ A2 cosh/3, 



(52) 



where ^3 = B/8 is to be determined, and A\ is given by eq. (1107p . One can fix ^3 using the require- 
ment that the potential must not diverge at (3 = in, if the spurious charge has been consistently 
removed: 



ITT 

The solution of the Newtonian potential which respects all boundary conditions is therefore 



(53) 



$(/3,r) = m 2 e T 



16vr 2 



(*/3 + 7T) 



1 



sinh (3 



+ 2 sinh j3 



+ C cosh (3 



(54) 



where C = ^2/772-2 plays the same role as the constant in the QED case. We see that the time- 
dependent potential contains the combination m2e T , as a consequence of the same combination 
occuring on the right-hand side of eq. (|42p . We can thus identify a static potential 



$(/3) = m 2 



16vr 2 



{i/3 + tt) 



1 



sinh j3 



+ 2 sinh (5 



+ C cosh (3 



(55) 



This is correct at all times if we consider mi to be time-dependent (777-2 (t) = "^2e T , where 777,2 is 
the static mass), which takes into account the fact that gravitational Wilson lines are not invariant 
under reparametrisations of p^. 



13 



Now consider adding a test particle of mass mi at location (3. If m 2 is at f3 = 0, then (3 = /3i2 is the 
geodesic separation between the masses in the AdS space. The total time-independent potential is 
then given by 



16vr 2 



sinh f3 



+ 2sinh/?i2 



.12 



+ C cosh /3i2 



(56) 



Note that the potential energy is defined only up to an arbitrary amount of the solution to the 
homogeneous equation^!. In the present case, this includes a contribution proportional to cosh {3, 
as seen explicitly in eq. (|46p . The constant C can be fixed, as in the QED case, by analytically 
continuing eq. (i56"j) to the case of one incoming and one outgoing particle, corresponding to the 
transformation of eq. (j50|) for the test particle of mass mi . Then 



E{(3 



12 J 



Wtt- 



12 



1 



sinh/3 



+ 2sinh/?i2 



12 



iC cosh f3 



12 



(57) 



This must vanish as (5 — ¥ 0, corresponding to the fact (in the original Minkowski space) that the 
anomalous dimension vanishes for a straight Wilson line contour with no cusp. This fixes 



C 



IK 

16^2' 



so that the potential energy in the original setup of two final state Wilson lines is 



E(j3 12 ) = m 1 m 2 



16vr 2 



(i£i2 + ?r) 



1 



sinh j3 



12 



+ 2 sinh P12 — i cosh /3\2 



Finally, the cusp anomalous dimension in Minkowski space must be given by 



(58) 



(59) 



r = iE = mim.2 



16vH 



(iir - /3 12 ) 



1 



sinh/3 



12 



+ 2 sinh /?i2 + cosh j3i 2 



(60) 



This agrees exactly with the field theory calculation of this quantity in Minkowski space, which we 
present here in appendix [Cl 

If many different masses are present, one finds a total anomalous dimension 



tot 



EK 
1fv 



16^ 2 1 



(in - fa) 



1 



sinh (3. 



+ 2 sinh (3ij I + cosh (3. 



i.i 



(61) 



corresponding to the total energy in the radial coordinate space. The sum is over all distinct pairs 
of outgoing particles, such that the cusp angle in each case is given by 



cosh f3i 



Pi ■ Pj 
miirij 



(62) 



Some further comments are in order regarding the anomalous dimension that we have derived. 
Firstly, note that eq. (|6ip is not the most convenient way to express this result, with a view to 



'in the QED case, for which Laplace's equation holds, this amounts to the presence of an arbitrary constant. 
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examining the light-like limit of — > oo. Instead, we may pull out an overall factor of cosh/3j 7 - in 



J i3 

each term to obtain the alternative form 

,2 



r 



K 



i<j 



16vr 2 



(wr - fiij) 2coth/% 



1 



sinh fiij cosh 



+ 1 



(63) 



where we have used the cusp angle definition of eq. (|109p . This will be useful in the following section. 



Secondly, there are potential conceptual issues regarding the gravitational cusp anomalous dimen- 
sion. This is defined in terms of the ultraviolet renormalisation properties of the vertex at which 
multiple Wilson lines meet. In calculating any particular diagram contributing to a Wilson line 
expectation value (and thus contributing to the gravitational soft function), one encounters addi- 
tional ultraviolet singularities relating in principle to the renormalisation of the masses rrtj, and the 
gravitational coupling constant At. However, gravity is non-renormalisable, such that the latter con- 
tributions are ill-defined, leading to additional operators at each order in perturbation theory that 
need to be included in Wilson line diagramJ^l. One may then worry that it is not possible to define 
a gravitational cusp anomalous dimension at all. However, this worry can be seen to be misguided 
for the following reasons: UV singularities of Wilson lines associated with renormalisation of the 
cusp correspond to IR singularities in scattering amplitudes, thus must be ultimately separable in 
a meaningful way from UV singularities associated with Lagrangian parameters. Also, as can be 
inferred from |50H52| . the gravitational cusp anomalous dimension is one-loop exact. Thus, the sit- 
uation in which one might have to worry about separating overlapping UV divergent contributions 
to the gravitional couplings and to the multieikonal vertex (which potentially occurs only beyond 
one loop order) never arises. It seems, then, that one-loop exactness of the gravitational anomalous 
dimension is crucially related to the fact that the IR limit of perturbative general relativity must 
be UV finite. 



4 Gravitational Wilson lines in the lightlike limit 

In the previous section, we have obtained the cusp anomalous dimension for perturbative GR via 
an energy calculation in radial coordinate space, similar to the QED / QCD analysis of [48J. To 
complete this analysis, it is instructive to examine the lightlike limit pf — > 0. As is well known, 
additional collinear singularities appear in this limit. In gravity, these singularities are present on 
a diagram by diagram basis, but cancel after summing over all diagrams [49U52] . The appearance 
of collinear singularities in the gauge theory case is associated in the radial coordinate space with 
a linear divergence of the (imaginary) energy of two static charges as (3 — > oo. A nice comparison 
was outlined in [48J between this phenomenon and linear confinement - the confinement in this case 
being that of e.g. quarks within jets. 

The aim of this section is to briefly revisit this analysis in the context of the gravity calculation 
carried out in the previous section, in order to complete the conceptual mapping between the gauge 
theory and gravity cases. We begin with the form of the anomalous dimension given in eq. (163p . in 
which an explicit factor of pi ■ pj has been pulled out in each term. This expression remains valid 
as mi — > 0, in which case fyj —> oo for all j. Let us now consider the case that particle i indeed 

15 For a recent discussion of the difficulty of considering running parameters in quantum gravity, see [70| . 
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becomes massless. From eq. (|63p . we may isolate all contributions involving particle i as 



Pi ■ Pj 



which as m,- 



(ft. 



16vr 2 
oo) becomes 



(iw - fa) 2coth/% 



sinh Pij cosh ftj 



+ 1 



(64) 



(65) 



Examining ft,- itself, this is 



ft 



cosh 



Pi • Pj 



log 



miiJij 

We may rewrite the logarithm on the right-hand side as 



Pi ■ Pj 
miirij 



log 



Pi ■ Pj 

rriimj 



log ( ) + log 



-Pi • Pj 

Qrrij 



Oil). 



log(^)+C)(l) 



(66) 



(67) 



where Q is an arbitrary momentum scale to keep the arguments dimensionless. Substituting this 
into eq. (j65j) gives 



r, 



8tt 2 



log 



Q 



Pi 



Using momentum conservation 



Pi = 0, 



the total contribution to the cusp anomalous dimension from particle i is 

„-2 



16vr 2 



log 



rrij 

Q 



Pi =o. 



(68) 



(69) 



(70) 



Thus, collinear singularities do not appear in the gravitational soft function. This is not, of course, 
a new result. The cancellation of soft collinear singularities has been known since |49| . and has 
been generalised to hard collinear singularities in [52J. The essential physical reason for this cancel- 
lation can also be obtained by direct analogy with abelian and non-abelian gauge theory: collinear 
singularities associated with a given particle depend on its squared charge. In QED this is q 2 , 
where q is the electromagnetic charge. In QCD, this is the quadratic Casimir invariant associated 
with the representation appropriate to the parton of interest. In gravity, the squared charge is 
the 4-momentum squared, which is zero if collinear singularities are to be present - which ends up 
removing them. 



Although the above argument is formulated for only one particle becoming massless, it generalises 
straightforwardly to cases involving more than one massless particle. Note that our reasoning does 
not tell us that massless particles do not contribute at all to the energy in radial coordinate space. 
Rather, the diverging term involving the geodesic separation of a given particle from one of its 
partners (as the former becomes massless) cancels after summing all contributions to the potential 
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energy, from all the other particles. There are still terms which are 0(/3°), which we neglected in the 
above analysis. It is in principle possible to calculate the total potential energy from interactions 
with the massless particle by giving the massless particle a small mass mi, as above, and setting 
mi — > at the end of the calculation after summing all contributions. 

How are we to square this with the fact that, for massive Wilson lines, we derived the energy in 
the radial coordinate space after taking the Newtonian limit? In this limit, it must be true that a 
massless particle contributes nothing to the energy of a collection of masses, as it does not gravitate. 
The resolution of this puzzle is that as m, — > 0, the Newtonian limit is no longer valid, as this relies 
on being able to define velocities which are much less than the speed of light, which is only possible 
for massive particles. Instead, as rrii — > 0, one must take into account special relativistic corrections. 

Having now completed our analysis of gravitational Wilson lines in AdS space, and their analogies 
with the QED case, it is amusing to note that both theories are in fact special cases of a general 
formulation, with a continuous relation between them. This is the subject of the following section. 



5 General formulation 

In the previous sections, we have reviewed the properties of Wilson lines in AdS space, and used 
an analogous analysis to [48J to examine the properties of the cusp anomalous dimension in per- 
turbative GR. In this section, we point out that one can formulate a general calculation for the 
cusp anomalous dimension, two special cases of which are QED and gravity. Furthermore, we will 
show that these cases are continuously related to each other. This perhaps adds an interesting ad- 
ditional way of thinking about the results of [48j and the present paper, which may be of further use. 



Consider the operator 



W n = exp 



iXp^p^ 2 ...p^ 



dsH. 



n (sp) 



(71) 



which is clearly defined for integer n. We may recognise this as a generalisation of the Wilson line 
operators of eqs. ([H l22j) . where p^ is the hard momentum of a particle which emits soft quanta 
of a spin-n gauge field H Ul ...^ n with coupling constant A 0. We have again parametrised the 
straight-line contour of the Wilson line according to eq. (|31h , which allows to rewrite eq. (171D as 



W n = exp 



exp 



/oo 
-oo 



(72) 



We see that the Wilson line phase is the product of a reparametrisation invariant quantity H Ul ^ an _ 1T 
(which has one radial-space index and n — 1 Minkowski-space indices) and a charge given by 



Xp 



.p^- 1 . The fact that the former is conformally invariant implies that 



H 



TT...T 



X^ 1 x^ 2 



X 



Hn-l 



MlA t 2---A t n-lT 



T H((3), 



(73) 



16 The reader may object, regarding how to interpret this operator for n > 3. We return to this point in what 
follows. 
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which is a generalisation of eq. (|30p . Let us now assume that the homogeneous equation of motion 
for H T ... T is 

DH T ... T = (74) 

(as is true for the QED and gravity cases). Implementing the behaviour of eq. (|73h and using the 
covariant D'Alambertian operator of eq. (|39p . one finds that the spatial part of H T ___ T satisfies 

[V 2 - (n 2 - 1)] E{P) = 0. (75) 

It is straightforward to verify that this equation reduces to eq. ([8]) and eq. (|44p for the QED (n = 1) 
and gravity (n = 2) cases respectively. The general solution of the general equation is given by 

where we have chosen to write this explicitly as the sum of two parts which have a definite parity 
under the transformation of eq. (|50p . One may check (using hyperbolic function identities) that 
this result indeed reproduces the QED and gravity results of eq. (f9"l I45|) . However, the solution of 
eq. (I76p (as a function at least) is well-defined for any n. In particular, we may consider n, divorced 
from the original context of eq. (|7ip . to be a continuous parameter that smoothly interpolates 
between the two solutions we obtained previously. 



In both of the cases considered so far, the solution of eq. (I76p did not have the right behaviour 
as j3 — > oo to correspond to the cusp anomalous dimension in the relevant field theory. This 
was rectified by modifying the current density in the inhomogeneous equation for the conformally 
invariant field component Hu 1 ... At „_ lT by a constant density. Introducing a delta function in the 
radial coordinate space, one may write the Wilson line phase from eq. (I72p as 



iXp 



/oo r 
dr / d*x5^(x)H^ n _ 1T {x), (77) 
-OO J 



such that the current density that sources the conformally invariant quantity n n _ 1T is given 
by 

^ = ^^ 1 ... Mn _ 1 T = ApM l _ _ . p Mn-l < j(3)( a .) j ( 78 ) 

from which one findj^l 

V-^n-ir = Ae- 2 > w . ..p^S^ix). (79) 

The appropriate generalisation of the constant charge density procedure is to modify this current 
to 



such that one finds 



* * * 1J -fJ.l...fJ, n -lT 

Xe^-^m 11 - 1 \5 {3) (x) + ETcosh^" 1 ) fi] . (81) 



7 Again we have absorbed geometric factors into the delta function, so that this is normalised according to eq. (|36|) 
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Thus, the constant charge density becomes a smooth distribution ~ cosh^ n ^ /3 distributed through- 
out space. As for eq. (I76D . one may continue n away from integer values. 

In the QED and gravity cases, the constant charge density procedure ensured that the energy asso- 
ciated with a pair of charges in AdS space diverged with an overall power of /3 at large separations 
(corresponding to collinear singularites in Minkowski space). This became somewhat non-trivial 
in the gravity example, in which it was crucial that the modification to j TT went like cosh/3. For 
integer values n > 3, we may note that this property does not generalise. One may verify that 

rv72 f 2 1 m ( (Ciff + C 2 ) cosh(n/3) \ _ 2Cm sinh(n/3) 
[V [H ijJ V sinh/3 J ~ sinh/3 ' 

r V 2 - (n 2 - 1)1 f ( C ^ + C ^ sinh (nft) \ = 2C 1 ncosh(n/3) 
L J \ sinh/3 / sinh/3 

Thus, a function constructed by modifying the solution of the homogeneous equation (eq. (|76p) to 
include an overall power of /3, is not consistent with a modified charge density ~ cosh^ n_1 ^ j3 in 
general (one may also show that the converse is true). Nevertheless, these results are correct in 
the QED and gravity cases. For QED, eq. ([82]) applies, and the resulting (constant) charge density 
picks out the solution of the homogeneous equation which has odd parity under the transformation 
of eq. (|50p . For gravity, eq. (|83|) produces the required cosh/3 charge density, which then picks out 
the even solution for H(/3). 

For n > 3, it is no longer true that the charge density in j T ... T arising from the constant charge 
prescription is such as to modify the energy by a linear term in /3. This presumably means that 
one cannot interpret the resulting energy as a cusp anomalous dimension of a Wilson line operator. 
However, this is not at all surprising, as the operator of eq. (|7ip ceases to be meaningful for n > 3. 
If eq. (|7ip is to be interpreted as the operator describing the emission of soft higher spin gauge 
bosons, then Lorentz invariance demands that the quantity 



be conserved, where the sum is over all external particles n, and g n is a constant which may depend 
on a given hard particle in general (see e.g. chapter 13 of [71]). The n = 1 and n = 2 cases 
correspond to electromagnetic charge and 4-momentum conservation respectively. However, for 
n > 3 no conserved quantity is possible if non-trivial scattering is to occur, which makes the above 
argument meaningless if n > 3 (for integer n). That this is seen directly in the general analysis, in 
terms of the charge densities not matching up, is interesting. 



In this section, we have considered a general formulation of the cusp anomalous dimension calcula- 
tion in Minkowski space, in which the QED and gravity theories emerge as special cases. Indeed, as 
discussed above, the general Wilson line operator of eq. (171 f) is only meaningful (for integer n) for 
n = 1,2. These are precisely the cases of QED and gravity respectively. Nevertheless, it is amusing 
to note that one may regard n also as a continuous parameter which smoothly interpolates between 
the one loop cusp anomalous dimensions of QED and GR. Our main motivation in pointing this 
out is in case this has any further interpretation or utility, in addition to novelty. 
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6 Conclusion 



In this paper, we have examined gravitational Wilson lines, representing the emission of soft gravi- 
tons from an eikonal emitter. We have examined in detail the anomalous dimension which controls 
the renormalisation of the vertex in a general correlator of such Wilson lines, with a view to gen- 
eralising the analysis of [48j which relates the calculation of such anomalous dimensions to static 
energies in Euclidean AdS space. 

There are a number of motivations for doing this. To start with, gravitational Wilson line operators 
and their renormalisation properties have only recently been studied. Given that the cusp anoma- 
lous dimension forms such a crucial object in gauge theories, any investigation of its properties 
in a gravitational context is interesting enough by itself. Our analysis also sheds further light on 
the radial coordinate picture for gauge theory Wilson lines developed in [IS], due to a number of 
subtleties in the gravitational case. For example, the procedure of adding a constant charge density 
to obtain the correct boundary conditions in the QED analysis implies a cosh (3 charge density in 
the gravity case, which is precisely such as to lead to the correct solution for the anomalous dimen- 
sion, and also to be consistent with the fact that masses, unlike electric charges, can only be positive. 

The similarities and differences between the gravity and QED / QCD radial coordinate pictures 
may be useful in providing further insight and intuition in current research at the boundary be- 
tween gauge and gravity theories. One example of this is the property of jet confinement in QCD, 
which has a particularly intuitive description in the radial coordinate picture, and which is absent 
in gravity. Another example where the radial coordinate picture might be useful is in thinking 
about the all-order structure of IR singularities. Radial coordinates were used in 08] to motivate 
a family of conformal gauges, in which the contributions from graphs involving multiple gluon ver- 
tices was much reduced (e.g. absent at two loop level). It seems likely that this is related to the 
conjectured dipole structure of soft singularities in QCD [30j[32j[33j[35]. In gravity there is less 
immediate motivation for considering such conformal gauges, due to the fact that the Wilson line 
operator is not conformally invariant, and also that the cusp anomalous dimension is known to be 
one loop exact. However, it may be that the simple structure of IR divergences in gravity (which 
is manifestly of dipole form), and the relationship between the radial coordinate pictures in both 
the gravity and gauge theory cases, can be used to gain further insight into the dipole conjecture. 

We also considered a general formulation of the cusp anomalous dimension calculation, in which 
a smooth parameter arises which smoothly interpolates between the QED and gravity cases. This 
analysis breaks down for n > 3 (where n is the number of momentum factors in the Wilson line oper- 
ator) , but this is entirely consistent with the fact that such operators are not physically meaningful 
due the fact that they lead to conserved higher-rank tensorial charges and thus trivial scattering. It 
is not immediately obvious what the applications of n considered as a continuous parameter might 
be, but one might hope that it may lead to further insight on the relationship between QED (or 
its non-Abelian brother, QCD) and gravity. 

In summary, the radial coordinate picture potentially offers novel new insights into both gauge and 
gravity theories, and the relationships between them. The results of this paper provide a useful 
addition in this regard. 
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(a) 



(b) 



Figure 2: (a) Wilson loop contour in Minkowski space used for the calculation of the Newtonian 
potential; (b) the relevant one loop diagram. 
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A Newton's law from the gravitational Wilson line 

In this appendix, we show how Newton's law of gravity can be calculated using the Wilson line 
operator of eq. (j22[) . First, we consider the Minkowski space contour C of figure [2{a). It is a 
textbook result in field theory (see e.g. [72]) that the static potential V(R) between two charges of 
separation R is given by 

V(R)= lim llog(W(C)>, (84) 

T— s-oo 1 

where (W(C)) is the vacuum expectation value of the Wilson loop along the contour C, and the 
latter is given in the present case by 



W(C) = exp 



jf dsp^p u h liu {sp >l ) 



(85) 



As the time T —> oo, the only contributing diagram at one loop order is that shown in figure [2](b), 
where the gravitons are emitted at positions 



x = spx, y = tp 2 , 
and the 4-momenta of the static masses are given by 

p 1 = (mi,0), P2 = (m 2 ,0) 



(86) 



(87) 
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One may evaluate this using the position space graviton propagator (in four dimensions) 

p. / _ \ _ _J PjJVjOfi 

V„v,ap{X V)- ijr2{x _ y y_ ie , 

as can be obtained by Fourier transforming the momentum space propagator of eq. (|19p . The 
diagram of figure [2]Jb) then gives a contribution 

W .0''f d s! Tlm ^ T ^^ + { ,% (89, 
\2J 4vr 2 J T/mi J (x-yy-ie 

where we have used the fact that the one-loop contribution to the Wilson loop expectation value 
is the same as the contribution to the exponent at this ordeiF^I. Using the definition of P^ v ,ap from 
eq. (fTUj) . one finds 

P m afi Pi Pi P2 P2 = 2(W • P2? ~ p\p\ = m\ml (90) 
Also transforming the integrals in eq (|89p to x° = sm\ and y° = srri2, one finds 

fti Tfl\TTli I o 



log W g = 2 / dx [) I dy 



16vr 2 J T J yv [-(x°-y°) 2 + R 2 -ie\ 
K 2 mini2T f T 1 

^OT-Tno 7yTT—> ( 91 ) 



~ 16vr 2 7_ T yu (y ) 2 -i? 2 + ie' 
where in the second line we have used the fact that we are taking T — > oo. The y° integral gives 



T ^ (y°) 2 -i? 2 +^~y-oo yu (yO_ jR + i e) (yO + ^_ ie) 



2 - ( ' (92) 



using Cauchy's theorem. Finally one finds 



= — iSfS- < 93 > 

Substituting this result into eq. (|84p and using the definition of k in terms of Newton's constant, 
k = ^16ttGn, gives 

V(R) = G ^ m \ (94 ) 
R 



(95) 



The force between the two charges is thus 

_ G N m 1 m 2 
R 2 ' 

which is Newton's law of gravity. 

18 In fact, the one-loop diagram enters the exponent to all orders, due to the gravitational analogue of Abelian 
exponentiation [49H52] . 
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Note that this calculation is closely related to an analysis carried out in [68], which studies non- 
perturbative quantum corrections to the Newtonian potential. To this aim, the authors define a 
Wilson line operator obtained by integrating over the worldline of a particle of mass m: 



W wi = exp 



-im I dt 
'c 



-9nu- 



dx^ dx v 
~dT~dT 



(96) 



where C is the worldline contour, and t the proper time along this path. This operator is equivalent, 
in the weak field limit, to the operator of eq. (|22p . as we now show. The weak field expansion of 
eq. (fT5j) allows us to rewrite eq. (f9"6"j) as 



wi = exp 



l J dtyj- (r]^ v + Kh fiu )pt 1 p u 



(97) 



where we have also taken the mass factor inside the square root and used p^ = mdx^ / dt. Expanding 
the square root to first order in the graviton field gives 



wl 



exp 



-im 



/ dt + i^ / dsh^p" 
Jc 1 Jc 



t 

m 



(98) 



The first term in the exponent is absorbed into the normalisation of the Wilson line operator. The 
second term is precisely that of eq. (|22p . That this must be the case follows from the fact that 
in the eikonal approximation, a particle emitting soft gravitons does not recoil, and thus follows 
its classical trajectory. Its action must then be given by its classical action, which is indeed the 
integral over the worldline as in eq. (|96p . This then fixes the form of the interaction between the 
eikonal particle and the (soft) graviton field. 



B Normalisation of the solution to Newton's equation 

In section^ we construct eq. (]42p for the Newtonian potential <3?, and show that the general solution 
of the homogeneous equation is given by eq. (|45p . In this appendix, we show how the constant A± 
can be related to the strength of the source term occuring on the right-hand side of eq. (J42]). First, 
note that rewriting this equation in terms of the covariant d'Alambertian operator gives 



1 



-.d u 



-9 



ni2e 



(99) 



Integrating over the complete 4- volume on both sides and using the covariant form of the divergence 
theorem |69|, one finds 



dn 



~im 2 e r (5 (3) (x), 



(100) 



where f2 is the 4-volume, and dS a the element of 3-surface area. The right-hand side gives 



-g m2 e 



'-) m 2 I dre 3r 



(101) 



where we have used the fact that y/—g ~ e 4r , and cancelled spacelike geometric factors in yj—g 
with similar factors in 5^(x), consistent with our definition of the delta function in eq. (|36p . We 



23 



T 





s 1 














S 2 









Figure 3: The 3-surface used to normalise the coefficient A\ in the general solution to the homo- 
geneous Newton equation. The horizontal plane is representative of the spacelike directions, where 
the radial distance in this plane represents /3. 



now consider the surface shown in figure El and consisting of a cylinder of radius f3o in the spacelike 
direction, with upper and lower surfaces at r = tq and r = respectively. Then eq. (jlOip becomes 



g) J dn^-gm 2 e^5^(x) = (|) m 2 J 



(102) 



Splitting up the surface into parts as labelled in figure [31 the surface integrals over S\ , S 2 and S3 
are given by 



r rPo 
/ dS^ y/^ggTdvQ = -47re 3ro / dp sinh 2 /3K(/3); 
JSi Jo 



__ 4-7T r 1 

dS„ ^g^d u $ = — sinh 2 /3 dp K((3) (e 3 ^ - 1); 

6 L J O 



IS 2 



S3 



dS„ d v <$> = 4vr / d/3 sinh 2 /3K(/3) 



(103) 
(104) 

(105) 



where we have used eq. (|43|) . Substituting explicitly the general solution of eq. (j4"5j) gives a total 
surface integral of 

r 4-7T 

dS„^g^d v $ = —A 1 (e 3T0 -l), (106) 



such that equating this with eq. (|1Q2[) gives 



l07T 



(107) 



as has been used in eq. ([4U|h 
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Figure 4: Diagram entering the calculation of the cusp anomalous dimension. 



C Field theory calculation of the gravitational cusp anomalous 
dimension 

In this appendix, we detail the calculation of the gravitational cusp anomalous dimension in 
Minkowski space, using a conventional field theory calculation. First, we consider the diagram 
shown in figure HJ consisting of a graviton exchange between the two contours on either side of the 
cusp. We consider the case that p\ and P2 are both outgoing, and correspond to different masses mi 
and m2- Thus, figure [4] is not a cusp in the traditional sense, but rather can be generally embedded 
into a graph where multiple contours intersect. Using the position space propagator of eq. 
this gives a contribution 



oo i 

dt- — . (108) 



o (sp 1 - tp 2 ) 



where the tensor is defined in eq. (fT9j) . In the (-,+,+,+) metric we are using for gravitational 

calculations, the cusp angle is given by 

cosh/?i2 = -ELJH (109) 

7711771,2 

where pf = —m 2 . Then eq. (|108p can be written, after transforming s — >■ s/mj, t — > t/m,2, 

F= _(_^j_ m^M r ds r dt i mo 

V 2) 4vr 2 mim 2 J J s 2 + t 2 - 2st cosh /3i2 ' v ; 

The integral over s and t is common to the QED case. However, we evaluate this here for 
completeness. First one sets s = Xt to give 



/ ds dt „ ,„ t . — = / -/ d\- ^4 — . (HI) 

Jo Jo 



**oc ) poo ^ poo 

s 2 + t 2 - 2st cosh /3i2 ~ Jo ~t Jo ""'l + A 2 - 2 cosh f3 12 ' 



The t integral contains an ultraviolet and infrared divergence, where the coefficient of the former 
gives the contribution to the cusp anomalous dimension at one loop order. At a given loop order, 
care must be taken to consistently separate the IR and UV singular parts, which becomes especially 
cumbersome when collinear singularities are present. At one loop, however, we may simply evaluate 
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the t integral by imposing UV and IR cutoffs (a similar procedure is used in the QED analysis 
of [IB])- That is, we may write 



00 dt 
T 



llH dt 



log 



Auv 



Air 
Auv 



Then one finds 



ds 



dt 



1 



s 2 + t 2 - 2stcosh(3 



12 



log 



A 



uv 



A 



1R 



dX 



1 



— cosh P12 



X 2 - sinh 2 # 



12 



(112) 

(113) 
cosh /?12- 

.„ (114) 

- cosh /3i 2 

Substituting eqs. (|113l 1114]) into eq. (|110p and taking the coefficient of log Auv > one finds a contri- 
bution to the cusp anomalous dimension given by 



where we have also completed the square in the A integral, and transformed A — > A 
Substituting A = sinh j3\2 coth u, the A integral can be carried out to givS 



dX 



1 



ITT 



12 



A 2 - sinh 2 p 



sinh j3 



Pi Pi PfMV, a p P2 P2 i7T ~ $vi 



16vr 2 



The kinematic factor is 



P2 P2 



m\m2 sinh j3\2 

'^{Pl-P2? 



(115) 



m\m2 



mi W2 
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mi rri2 [2 cosh 2 f}\i — l] 
mi?7i2 [l + 2 sinh 2 fiyi\ 



Finally, one has 



167T' 



: mi m 2 (i7r - /3i 2 ) 



1 



sinh /3 



+ 2sinh/3i 2 



12 



(116) 



(117) 



This is not the whole story. One must also add self-energy diagrams associated with each external 
line. Rather than calculate these directly, one can surmise their contribution as follows. Firstly, the 
effect of each self-energy diagram can only depend upon the quantum numbers of a single parton 
leg, so that the sum over all self-energy contributions has the form 

E^ 2 = E^'^ 



8=1 



(^2 Pi] -l^Pi-Pj 

\i=l / j>i 



(118) 



where C is a constant independent of the parton index i, and we have rewritten the momentum 
dependence on the right-hand side. Using momentum conservation (eq. (|69p ). one may rewrite 
eq. ([iTHI) as 



E^ 2 



-2C 



i=l 



i j>i 



(119) 



19 Care is needed with the minus sign in this equation, where the in results from correctly implementing the it 
prescription in the graviton propagator. Alternatively, one may carry out the calculation for one momentum incoming 
and one outgoing, before analytically continuing /3i2 = in — P12 ■ 
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and one sees that each pair of external lines is associated with the contribution 

—2Cpi ■ pj = 2C mi rrij ( — — ^ 

\ rrii rrij 

= 2C rrii rrij cosh , (120) 

where we have used the cusp angle definition of eq. (j62|) . Adding this to eq. (|117p , one may fix the 
constant C by requiring that the cusp anomalous dimension vanishes at /3i2 = ivr, corresponding 
to a straight line Wilson contour with no cusp. The complete result is then 



K 2 



( sha |^ 2 +2suih/j 12 j + cosh( )' l2 ) 



(121) 



16vr 2 

Note that in the large j3 limit, this has the behaviour 

r(/3) ~ PeP, (122) 

as discussed in section [3J 
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